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A new agent-based, bounded-confidence model for discrete one-dimensional opinion dy- 
namics is presented. The agents interact if their opinions do not differ more than a toler- 
ance parameter. In pairwise interactions, one of the pair, randomly selected, converts to 
the opinion of the other. The model can be used to simulate cases where no compromise 
is possible, such as exclusive choices or competing brands. The fully-mixed case with 
^ 1 maximum tolerance is equivalent to the Gambler's Ruin problem. A fully-mixed system 

' always ends up in an absorbing state, which can have one or more surviving opinions. 

' An upper bound for the final number of opinions is given. The distribution of absorp- 

, tion times fits the generalized extreme value distribution. The diffusion coefficient of an 

■ opinion increases linearly with the number of opinions within the tolerance parameter. 

• ' A general master equation and specific Markov matrices arc given. The software code 

developed for this study is provided as a supplement. 

Keywords: Opinion dynamics; diffusion; bounded confidence; agent-based modeling; cel- 
lular automata; non-compromise model. 



1. Introduction 



' Opinion dynamics, as the term is used in the physics community, studies the evo- 

lution of systems composed of a large number of individuals characterized by a 
state (the "opinion") and modify their state by repeated interactions whose rules 
are based on social observations rather than physical laws. The goal of opinion 
dynamics is to develop models that capture the essential behavior of social systems. 

Several models for opinion dynamics are developed in the last decade 6J • Some 
of them are agent-based, containing a large number of simulated persons, called 
agents. Depending on the details of the dynamics, they may all converge to a single 
opinion (consensus), to two separate opinions (polarization) or to several opinions 
(fragmentation). Usual questions include whether agents always (or at all) reach a 
consensus, and under which conditions; how long it takes to reach an equilibrium 
state; what the effect of the topology of the interaction network is, etc. 

In models where agents are chosen randomly in pairs, the updated opinion of 
each agent is a function of the interacting pair only. Examples include the Axel- 
rod model 0], the voter models [H [55], the CODA model [HI [H], the Deffuant- 
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Weisbuch model [9l[30], and their variants. The differences and similarities between 
these models and the present model are discussed in section [TT2] 

In the so-called "bounded confidence" models, two agents interact only if the 
difference between their opinions is less than a given threshold value, called "uncer- 
tainty", "tolerance", or "confidence bound". This constraint reflects the fact that 
people change their opinions only to a limited degree. Each person has a certain 
tolerance interval, and ideas that are outside this interval are deemed too radical; 
they are not even considered or discussed. Bounded confidence can be applied to 
either type of models discussed above, but it makes sense only if there are three or 
more opinions. 

This paper is organized as follows: Section 11.11 introduces the model, and ex- 
plains the rationale behind its assumptions. Section 11.21 compares the model with 
existing models that use pairwise interactions. Section [2] describes how agents move 
in opinion space with or without bias, and outlines the general properties of the 
absorbing (final) states. Section [3] analyzes fully-mixed populations using Markov 
matrices and agent-based simulations. Section |4] briefly discusses some features of 
the model under general network topologies. 

The C and Matlab codes of the programs used in this study are provided as a 
supplement to this paper. Those who wish to replicate the presented results can 
download the codes from the ancillary files link in ArXiv, or from: 

https : / /sites . google . com/site/mkaaiiozturk/progrELins/opdyncode .zip 

1.1. The model of interaction 

The model comprises N agents, each carrying an integer opinion in the range [1, Q]. 
One agent and one of its neighbors in the network are chosen randomly, with opin- 
ions q and g', respectively. They interact only if 1 < \q — q'\ < rf, that is, if the 
opinions differ by no more than the confidence bound d, an integer between 1 and 
Q — 1. If this condition holds, we call the pair of agents "compatible". 

In this study, all agents in a given population have the same value for d. Also a 
completely connected interaction network is assumed (except for the last section), 
meaning that any two agents can be paired up. 

At every time step of a simulation one pair of compatible agents is chosen ran- 
domly, so that at every step there is exactly one opinion change. Thus, "time" 
is defined as the number of conversions, not as the number of random meetings 
between agents. 

Upon interaction, a random choice is made so that agents both adopt opinion q 
with probability p or opinion g' with probability 1 — p. In the symmetric (unbiased) 
case p is set to 0.5. 

A second random choice may look superfluous, as agents are already chosen 
randomly. However, bias can be introduced by making the probability p a function 
of the properties of agents. For such second random choice is necessary at 

the interaction stage. 
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Dynamics of discrete opinions without compromise 3 

Some opinion dynamics models (e.g. Deffuant-Weisbuch j9]) prescribe compro- 
mise: Opinions of both agents are updated so that their opinions come closer to each 
other. However, such a compromise is not always possible. When a couple disagrees 
in what movie to see, they do not watch half of each. Supporters of two electoral 
candidates do not agree to vote for a third candidate. When a Mac user and a PC 
user interact, they can rarely agree to construct a hybrid computer. The presented 
model applies to situations where a compromise is not possible or reasonable. 

The lack of compromise implies discreteness of opinions because the initial set 
of opinions in the population is closed under the described interaction. However, 
this discreteness is more than a technical necessity; it is a fundamental requirement 
dictated by the limits of human cognition. 

In this study we deliberately use a small number of opinions. Many instances 
in life present us with just a few discrete choices which are set beyond our control. 
We vote for one of the existing political parties, or we purchase one of the available 
brands of toothpaste. Even if do not like any of these choices, we usually cannot 
construct the choice that we prefer. Furthermore, having too many choices is not 
necessarily better: Consumers experience "information overload" when faced with 
more than ten options for a specific item, resulting in confusion about choices and 
dissatisfaction with their final choice |15| . 

It can be argued that neither continuous opinion values nor a large number of 
discrete choices are consistent with the real decision-making mechanisms of the hu- 
man mind. Observations about the throughput of the human mind, reviewed by 
G.A. Miller [5D], suggest a realistic upper limit for the number of choices. In his 
classic paper. Miller reviews experiments in which people are asked to differentiate 
between various one-dimensional sensory stimuli, such as sound pitch, sound inten- 
sity, shape sizes, colors, brightness, and vibration intensity on skin. Miller concludes 
that the "channel capacity" of the human mind is limited so that we can distin- 
guish only about 5 to 9 different categories of a one-dimensional stimulus, hence 
the famous "seven plus or minus two rule" . 

For the same reason, questionnaires used in social research involve interval scales, 
such as the Likert scale, consisting of only five or seven (or occasionally nine) lev- 
els of responses [T9|. Studies indicate that respondents cannot distinguish further 
subdivisions of attitudes and opinions [121 Hj ■ These observations suggest that an 
opinion dynamics model must not have too many opinions, lest it disagree with 
human cognitive limits. 

In a similar reasoning, Urbig 12 7j points out that an attitude and its communi- 
cation are different things, and criticizes models that are "based on the assumption 
that individuals can communicate the difference between an attitude of 0.5555 and 
0.5556.'''' In Urbig's sophisticated model, agents have continuous opinions but when 
they interact only discrete values are communicated. 
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1.2. Comparison with other models 

In this section the model described above is compared to existing similar models. 
Most of these models have been extended in different levels of sophistication, e.g. 
using heterogeneous agents or complicated network topologies. Here we consider 
only their fundamental features. 

In Axelrod's model for the dissemination of culture [3], agents' opinions are 
vectors composed of integers. At each step one agent (called the active site) and 
one of its neighbors are selected randomly. The two interact with a probability 
proportional to the number of opinions they agree on. If interaction takes place, the 
active agent copies one of the differing opinions of the neighbor. 

The only common feature between the Axelrod model and the current model is 
that both are non-compromising. One agent adopts the opinion of the other as it is, 
instead of both of them converging to each other. The Axelrod model does not use 
bounded confidence; instead, the proximity of opinions is measured by the number of 
identical opinions. Also, the Axelrod model makes sense only with vector opinions. 
A full comparison between the models requires extending the presented model so 
that agents hold a vector of opinions. This is a straightforward generalization, but 
out of the scope of this study. 

The model by Laguna et al. |14) uses vector opinions with binary values in 
each component. It adds bounded confidence to Axelrod's model by making two 
agents interact only if the Hamming distance between them (the number of different 
components) is less than a threshold. 

The classical Voter Model [6] assigns a scalar binary opinion to each agent. One 
agent is chosen randomly and it adopts the opinion of a random neighbor. This 
is a very simplified version of the Axelrod model. As there are only two opinions, 
bounded confidence does not apply and any agent can interact with any other. The 
Constrained Voter Model 28 extends it by introducing an intermediate opinion. 
Agents can be leftists, centrists or rightists. Leftists and rightists do not interact 
with each other, but with centrists only. These models are special cases of this 
paper's model: The classical Voter Model is obtained by setting (5 = 2, the Con- 
strained Voter Model is obtained by setting Q = 3, d = 1. 

Mobilia [21] introduces a control parameter a to the Constrained Voter model 
so that when a > the absorbing state is more likely to be composed of extremists, 
and when a < it is more likely to be centrist. This can be implemented in the 
current model by making the probability of conversion p a function of interacting 
opinions. 

The Deffuant-Weisbuch Model [HI 130] j in its basic form, uses continuous opinions 
in the range [0, 1]. At each time step, one agent is chosen randomly, and its opinion 
value X is compared to a randomly chosen neighbor's opinion x' . If the absolute 
value of the difference between opinions is less than the tolerance, the opinions of 
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agents are updated according to: 



X <— X + fi{x' — x) 

a;' •(— x' + — x') 



(1) 



where the convergence parameter fi G (0, 0.5] specifies how much the opinions ap- 
proach each other. The model imposes some compromise between agents. After 
interaction, each agent is more predisposed to its partner's opinion. When fi = 0.5, 
both agents agree on a common opinion, the mean of previous values. 

Discretized versions of the model round the opinions to the nearest integer 
(e.g. [IS]). In case of binary opinions, one agent adopts the opinion of the other 
with probability ^ [9l [26]. Alternatively, if opinions differ only by one, one agent 
can be made to adopt the opinion of the other agent with probability 0.5 [5]. This 
rule would make the model equivalent to this paper's model with d = 1. 

Although there is not a simple equivalence between the Deffuant-Weisbuch 
model and the present model, they have some similarities. In both models consen- 
sus is more likely as the tolerance parameter increases. Simulations [9l [11] indicate 
that for continuous opinions, the population almost always ends up in consensus 
when the tolerance parameter is larger than 0.5, or in discrete terms, larger than 
Q/2. In the present model, consensus is the only possible outcome if and only if the 
tolerance parameter d is Q — 1, its largest possible value. 

Agents in the Continuous Opinions - Discrete Actions (CODA) Model [16] carry 
only binary (or ternary [17] ) values of opinions, and an agent converts to the rival 
opinion with a probability that changes in time. This probability is equal to the 
degree of belief of the agent that the other opinion is true. At each step the agent 
observes the opinion of its partner, and updates the belief using Bayes' theorem. 
Effectively, agents have a memory of past interactions. The CODA model is similar 
to Urbig's model [37] in its observance of verbalization limits. In this paper's model 
currently there is not an individual conversion probability associated with each 
agent, so it is not as general as the CODA model. 
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2. Dynamical features 

We define the popularity Uq of opinion q as the number of agents with opinion q. 
Figure [1] shows a typical time series of popularities of five opinions. The opinions 
are initially distributed uniformly over 100 agents. The confidence bound d is equal 
to one, and there is no bias in the interactions (both agents in a chosen pair are 
equally likely to be converted). In this particular smulation we see that opinion 4 
loses all of its followers by step 300. After that time agents with opinion 5 cannot 
interact with any other agent, so the number of agents with opinion 5 does not 
change anymore. Similarly, agents with opinions 1, 2 and 3 keep interacting until 
the number of followers of opinion 2 drops to zero at step 1017. After that, no more 
interactions are possible as the surviving opinions 1, 3 and 5 differ by two, which is 
more than the confidence bound d = \. 

2.1. Agents' motion 

The agents move randomly in the opinion space as they interact with each other. 
Even if an agent starts with an extreme opinion, in time it can move to farther, 
initially incompatible opinions, as long as there are other compatible agents with 
whom it can interact and use as a stepping stone. 

An opinion that loses all of its followers is forever lost (unless the model in- 
volves some random noise that repopulates it). Absence of an opinion may block 
the migration of agents in the opinion space. For example, if there are no agents 
with opinion 3 and if the tolerance parameter d is equal to one, there will be no 
crossover between opinions 2 and 4 because such an interchange would require in- 
teraction with opinion 3 first. If d = 2 such a crossover is possible as agents with 
opinions 2 and 4 can directly interact. 

2.2. Relation to the Gambler's Ruin problem 

With two opinions, the dynamics is equivalent to the classical Gambler's Ruin prob- 
lem [lOj . where each player gains or loses one unit with equal probability. Here, play- 
ers are opinions that gain or lose followers. Even if one starts with several opinions, 
toward the end of the run, two opinions may be left within each other's confidence 
bound, reducing the system to the classical problem. 

The multiple player Gambler's Ruin is described by Ross [23] as follows: Q 
players (opinions) begin with different amounts of capital (popularity) riq {q = 
1 . . .Q) at their disposal. At each step two players are chosen to play one round 
and each player is equally likely to win one unit. Any player whose capital drops to 
zero is eliminated. The run continues until all units are accumulated by one player. 
So, this game is equivalent to this paper's model with unbiased interactions and 
d = Q-l. 

Under these conditions, regardless of how pairs of players (opinions) are chosen, 
Ross proves that: 
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Dynamics of discrete opinions without compromise 7 

(1) The expected number of interactions is finite, and given by (^N^ — "-g) /2- 

(2) The probabihty that opinion q is the ultimate winner (consensus opinion) is 
njN. 

(3) The expectation value of the number of interactions involving opinion q is 

nq{N - Uq). 

(4) The expectation value of the number of interactions involving only opinions q 
and q' is UqUqi . 

The related "iV- Tower problem" considers the game ended when only one of 
the players loses all of its capital. Bruss et al. [5] give the probability distribution, 
mean and variance of the absorbing time for the case of three towers (opinions). 
The authors argue that simple expressions may not exist for these quantities if more 
than three opinions are present. 

2.3. Absorbing states 

After repeated interactions the system eventually reaches an absorbing state where 
no further change is possible. This can happen in two ways: 

(1) No agent is compatible with any of its neighbors. This case happens when the 
network of agents has low connectivity (a few connections per agent). It is 
discussed in more detail in Section |4l 

(2) An opinion q has nonzero popularity and all compatible opinions q — d,...,q + d 
have zero popularity, so no interactions are possible. When the agents' network 
is completely connected, this is the necessary and sufficient condition for sta- 
tionary state. 

If there is a bias in the interaction that favors less-popular opinions, or if there 
is random noise that repopulates empty populations, the system may not be able 
to arrive at an absorbing state. 

We call an opinion stationary if all opinions compatible with it have zero pop- 
ularity. In the unbiased case and with some biases (see subsection 12. 4|) the system 
almost surely ends up in an absorbing state composed of several stationary opinions. 
Each stationary opinion is separated by at least o? -I- 1 empty opinions. Because of 
the stochastic nature of the dynamics, neither the stationary opinions themselves 
nor their number can be predicted in advance. The same initial conditions may 
lead to a different set of end states in each simulation. However, we can determine 
the maximum number m of stationary opinions by considering the tightest possible 
arrangement of stationary states. If Q is the total number of opinions and d is the 
confidence bound, starting from opinion one, the tightest stationary arrangement 
is obtained with steps of c? -I- 1 : 

1, l + (rf+l), l + 2(rf+l), ■•• l + (m-l)(d+l) <Q. (2) 

Alternatively, one can start at opinion Q and make another tight arrangement by 
going backwards with steps of c? -I- 1. Either case has m states, where m is the 
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maximum number of stationary opinions, given by 



1 



1, 



(3) 



where [-J denotes the integer part. In particular \i d = Q — 1, only one stationary 
opinion (consensus) may exist. For d < Q — 1, there may be one or more stationary 
opinions. 



2.4. Biased interactions 

In order to make the model more realistic we can introduce a bias in the interaction. 
This can be achieved by making the conversion probability p as a function of agent 
properties. 

When agents are identical, p can depend on interacting opinions g, q' . With 
proper choice of parameters, such a bias can reproduce Mobilia's |21J version of the 
constrained voter model. 

Alternatively, the bias could depend on the popularities n, n' of interacting 
opinions q, q' , respectively. For example, if an opinion is shared by more people, 
it could gains new followers more easily. This version of Matthew effect can be 
modeled in many ways, including: 

(1) Pairwise majority bias: If n > n', q is favored with probability 

p = 0.5 + 0.5\n-n'\/{n + n'). (4) 

(2) Global majority bias: If n > n', q is favored with probability 

p^0.5 + 0.5\n-n\/N (5) 

where N is the total number of agents. 

The pairwise majority bias considers only the relative difference between the two 
interacting opinions. The global majority bias puts this difference in a system-wide 
perspective. With the latter, the unfavored agent is not much impressed if the other 
agent's opinion is not very popular over the entire population. The difference be- 
tween these two biases is significant only when there are many competing opinions. 
Any bias to majority will accelerate the convergence to an absorbing state. 

Strictly speaking, as the popularity is a global variable, its use is not compatible 
with mindless-agent based simulation. Still, in a completely connected society where 
everybody knows everybody, it is not unreasonable that every agent has an idea 
about the popularities of opinions. 

Figure [2] compares typical time series under different bias conditions. Under 
pairwise majority bias, the time series is much smoother and reaches the end state 
rapidly. 

One can also set up a bias to minority opinions. For example, if the popularity of 
an opinion is too high, agents may be set up to "seek novelty" so that interactions 
are biased toward less popular opinions. In such models stationary opinions may 
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Fig. 2. A comparison of time evolution of popularities under no bias, pairwise majority bias and 
global majority bias. 5000 agents, 20 opinions, confidence bound 1. Only some opinions arc shown. 

not exist, but a dynamical equilibrium may be possible. Such biases are not studied 
in this paper. 

3. Fully mixed system 

In a "fully mixed" system, all agents have the same features and pairwise inter- 
actions between agents do not depend on individual properties. Furthermore, any 
agent can interact with all other agents as long as the bounded confidence con- 
dition is satisfied. A fully mixed system is a model for a small community where 
individuals are free to move around and discuss. 

3.1. Markov chain analysis 

Let Xt be the state of the system at step t, and let Pr{Xt = (ni, n2, . . . , ng)} be 
the probability that at time t opinion 1 has rii followers, opinion 2 has 712 followers, 
etc. By the law of total probability we have: 

Pr{Xt+i ^s} = Y, MXt+i = s\Xt = s'}Pv{Xt = s'} (6) 

s' 

where the sum is over all states. At each time step a pair of agents with compatible 
opinions q and q' {q ^ q') are selected. There are only two possible transitions: 
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nq ^ Hg — 1, Uqi Uqi + 1, Or ,nq Uq + 1, Uqi — )■ Uqi — 1. For both of these 
transitions, the transition rates are given by: 

where p = 0.5 for the unbiased interaction, and 

Q-l min(fe+rf,Q) 

fe=i fe'=fe+i 

is the number of all compatible pairs of agents in a given configuration rij. Note 
that M{{ni}) — for absorbing states. 

Let s be the state that the popularity of opinion i is n^. Let q, q' be two given 
opinions and let s' be the state where the popularity of g is — 1, the popularity 
of q' is Uqi +1, and the remaining popularities are the same as in s. Then the 
stochastic dynamics of opinion popularities is described by the following coupled 
linear system: 

Pr{Xt+i = s} = 

Q Tam{q+d,Q) 

^ P Mis') ~ l)^''^^* = ^ > 

g=l g'=max(l,(j— d) 

+ <5M(«),oPr{Xt = s) (9) 

The sum is over all ordered compatible opinion pairs ((7, g') where (7' ^ q. The step 
function 0(ng — 1) ensures that opinion q has at least one follower. The last term 
sets the absorbing state probabilities using the Kronecker delta symbol. 
With N agents and Q opinions, the number of possible states is 

+ (10) 

which is also the number of equations. With a proper enumeration of states, this 
linear system can be converted to a matrix- vector equation. The matrix is a Markov 
matrix, so all of its eigenvalues are less than or equal to 1. The eigenvectors corre- 
sponding to unit eigenvalues are absorbing states. 

The eigensystem of the Markov matrix can be determined mimerically, yielding 
the relative probability of each absorbing state, as well as transient properties like 
absorption times. However, as the matrix size grows as 0{N'^~^), this approach is 
feasible for small systems only. A direct agent-based simulation is needed for larger 
systems. 

The probabilities are expressed so that there is exactly one conversion between 
opinions, in accordance with the agent-based simulations. It is possible to modify the 
model so that at each time step any two agents are chosen randomly, compatible or 
not. In that case the system will have the same absorbing states, but its convergence 
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will be much slower. For this modification it is sufficient to set M{s') = N(N — 1) 
in Eq.Q, and replacing the last term of the same equation with 

3.2. Three opinions 

Here we apply Eq. ([9|) to the case of three opinions in a population of N agents with 
unbiased interactions. 

Let the confidence bound d be 1. By inspection it is seen that there are + 2 
absorbing states: (0, N, 0), and (i, 0,N — i) where i = . . . N . Compatible opinion 
pairs are (1,2), (2,1), (2,3) and (3,2), and M{x,y,z) = xy + yz. Then, Eq.® takes 
the form: 

Pr{Xt+i = {ni,n2,n3)} = ^ 
+ 



1 (ni 
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+ l)e(ni - 


1) 
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"-2 - 1,^3 + 1)} 


Pr{Xt = 


("1 


"-2 + l,"-3 - 1)} 



+ ^M(ni,«2,n3),oPr{Xt = (ni,n2,»^3)} (12) 

As 713 = — ni — 712, a state can be specified with ni and 7^2 alone. We list the 
states such that 77,2 varies from to TV — ni, and ni varies from to N, and we 
define i(7ii, 77,2, 773) to be the position of state (711, 712, 773) on this list (77,3 is kept for 
notational convenience). The list goes as follows: 

i(0,0,iV) = 1 
i(0,l, A^- 1) = 2 



i(0,iV,0) + 1 

i(l,0, A^- 1) = A^ + 2 

t{N,0,0) = (A^ + 2)(A^ + 3)/2 (13) 
Then, the index function must have the following form: 

i(77i, 772, n^) ^ {N ~ 2)77.1 - ^rii{ni + 1) + 772 + 1 (14) 

Now we define a column vector p* = [p*] such that = Pi'lATf = 
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(ni, n2, n^)}- Then the Markov matrix A = [a.y] is given by 

^^(ni.na.na) ,j — 9 _ 1 J_ n J_ -X- 1 '>n„ -l'"2 + l:n3),J 

'^j(rii + l,n2-l,n3),J 

^^2(711 ,712 — l,n3 + l),J 



+ ^i(ni,n2,n3) .i'^Af (ni,n2,n3),0 (l^) 

SO that p* = Ap*~^ 

Let Ai,A2,... be the eigenvalues of arranged in decreasing order, and let 
vi, V2, . . . be the associated eigenvectors. As the system has A^ + 2 absorbing states, 
it holds that Ai = . . . = AAr_|_2 = 1, and all other eigenvalues are smaller than unity. 

Matrix A is diagonalizable, so the equation can be written as: 

p* = VDV-^^^-^ = FDV^^p" (16) 

where V is the matrix whose columns are the eigenvectors of A, I? is the diagonal 
matrix whose entries are the eigenvalues of A, and p° is the initial probability vector 
of states. The vector l/~-'^p" = c gives the initial vector in the eigenvector basis. 
Then it holds that 

p* = V^L>*c 

= A*iCiVi+A^C2V2 + ... (17) 

where Cj is the i-th component of c. As t — >■ oo, all terms with A.; < 1 will vanish, 
and the steady state probability vector will be a linear combination of absorbing 
state eigenvectors: 

P°° = CiVi + C2V2 + . . . + CAr+2VAr+2 (18) 

Therefore, Ci is the probability that the system ends up in the absorbing state 
corresponding to v^. Numerical results for these absorbing state probabilities are 
given section [231 

If the confidence bound d is 2, compatible opinion pairs are (1,2), (1,3), (2,1), 
(2,3), (3,1) and (3,2), and M(x, y, z) = xy + yz + zx. Then the Markov matrix is 
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given by: 

0'i{ni,n2,n3),j = 7^ ^\(^^ j_ ^\ ^ 1^^ _L 1\^^ A_ ^-^^^^Oiin^-Y^n^ + l^n^),] 

" <5i(rii +1,712 — 1, "3), i 
^i{n\ ,n2 — 1,713 + 1), J 
^i(n\ ,n2 + l,n3-l),j 

'^i(rii+l,n2,n3-l),i 
<5i(rii-l,n2,n3 + l),j 

~^ ^2(ni ,n2 ,n3),j ^M(ni ,n2 ,713) ,0 (-^'^) 

which can be analyzed as described above. 

Using this method, one can set up expUcit Markov matrices for arbitrary Q 
and d. However, the explicit form of the index function that maps the 

linear index to states is too difficult to determine for general Q and N . Instead of a 
functional form, one can use a lookup table for this map when setting up the matrix 
in software. 

Matlab code that evaluates the absorbing state probabilities for Q = 3 and 
(5 = 4 using Markov matrices is given in the supplement to this paper. 
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+ 1) 



3.3. Absorbing state probabilities: Consensus or polarization? 

The system is said to be in consensus if all agents finally end up in the same opinion, 
polarized if there are only two remaining opinions (necessarily incompatible with 
each other), and fragmented otherwise. 

For Q — 2, only consensus is possible. If Q = 3 and d ~ 1, consensus or po- 
larization is possible. Fragmentation (three or more stationary opinions) is possible 
only for Q > 5. 

In general, there can be at most [{Q — l)/((i+ 1)J + 1 opinions in the absorbing 
state. If d = Q — 1 there can be only consensus. If (Q — 3)/2 < d < Q — 1 holds, 
the system can reach consensus or polarization, but not fragmentation. 

Using the Markov matrix, the probability of an absorbing state can be deter- 
mined as described in Section [321 For Q = 3 and d — 1, possible absorbing states 
are (0, N, 0) and {i, 0,N — i), where i = . . . N. Using a uniform initial condition 
(rti = 712 = n3 = -^/3), the matrix is numerically analyzed for N = 15,30,60, and 
99. In each case, it is found that the probability of the absorbing state (0, N, 0) is 
1/3, within numerical precision, independent of TV. 

Figure [3] shows the probability of the absorbing state {i, 0,N — i) as a function 
of z, for various N , for unbiased and biased interactions. Bias is taken to be pairwise 
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Fig. 3. The probability the absorbing state (i, f),N — i) for Q = 3, o! = 1. Left: Without bias. 
The curve for A'^ = 300 is obtained by direct agent-based simulation, the rest by Markov analysis. 
Right; With bias to pairwise majority (Eq. |4j; all curves obtained by direct simulation. 



majority bias, as given in 

In the unbiased case (left panel in Fig. ^ the consensus states (0,0, A^) and 
(A^, 0, 0) are more likely than polarized states. The probability of consensus states 
is about 0.065, and varying only slowly with N . 

The case of = 300 could not be analyzed with the Markov approach as the 
matrix is too big. This curve is obtained by using direct agent-based simulation, av- 
eraging over 10"* runs. In the biased case all calculations are agent-based, averaging 
over 10^ runs. 

In all cases the curves have a peak in the middle, indicating that an equal division 
{N/2, 0, N/2) is the most likely polarization. In the unbiased interaction the curve 
becomes flatter as N increases, suggesting that all polarization states are equally 
likely as iV — ^ oo. In the biased interaction, however, the central peak persists, 
therefore the near-equal polarizations are dominant. 

Table [1] lists the probabilities of finding one opinion (consensus) , two opinions 
(polarization), and three and four opinions (fragmentation) in the absorbing state, 
for Q up to 7. The probabilities are found by agent-based simulation, averaging over 
lO'"' runs. The system is initialized with equal popularities for all opinions. Biased 
simulations use pairwise majority bias. The number of agents is either 100, or the 
largest number smaller than 100 divisible by Q. For any Q, runs with d = Q — 1 
always result with consensus, so they are omitted in the table. 

The results are not sensitive to the number of agents. Larger populations (up 
to 1000) change the values only by about 0.01, usually less. 

The constrained voter model as described by Vazquez and Redner 28] is equiv- 
alent to the case Q = 3, d = 1, unbiased interactions. The authors give the prob- 
abilities of polarized states and consensus states as functions of initial densities. 
Applying these results to the uniform initial distribution, one finds that the prob- 
ability of a polarized state is 0.5377, probability of consensus on either extreme is 
0.0645, and the probability of consensus on the middle opinion is 1/3. These values 
agree with the results given above and in Table [TJ 
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Table 1. Probabilities of the number of opinions in the absorbing state with N ^ 100 agents, 
averaging over 10^ runs. 



no bias 


1 


2 


3 


4 


0.4634 


0.5366 






0.1255 


0.8745 






0.6735 


0.3265 






0.0236 


0.6770 


0.2995 




0.4042 


0.5958 






0.7831 


0.2169 






0.0029 


0.3309 


0.6662 




0.1897 


0.8103 






0.5680 


0.4320 






0.8425 


0.1575 






0.0002 


0.1149 


0.7151 


0.1698 


0.0839 


0.7878 


0.1283 




0.3832 


0.6168 






0.6731 


0.3269 






0.8842 


0.1158 







with bias 


1 


2 


3 


4 


0.4644 


0.5356 






0.0079 


0.9921 






0.6142 


0.3858 






0.0001 


0.6294 


0.3704 




0.3134 


0.6866 






0.6967 


0.3033 






< 10-5 


0.1854 


0.8146 




0.0216 


0.9784 






0.4649 


0.5351 






0.7529 


0.2471 






< 10-5 


0.0116 


0.7474 


0.2410 


0.0032 


0.7714 


0.2254 




0.2532 


0.7468 






0.5615 


0.4385 






0.7913 


0.2087 







The chance of reaching a consensus varies greatly with the confidence bound 
d. When the confidence bound is smah, the population is more hkely to end up in 
polarization, or even fragmentation, instead of reaching a consensus. This is a com- 
mon feature of bounded-confidence models: It is harder to reach a consensus when 
individuals do not have a high tolerance for different opinions. Narrow-mindedness 
splits society into several factions that do not speak to each other. 

Still, the population is usually not as much fragmented as it can possibly be. 
The entries for Q = 5 and Q = 7 with d = 1 show that the probability of the most 
fragmented state is smaller than a less-fragmented one. 

When interactions involve bias to pairwise majority, the same properties hold, 
but the system tends to be slightly more polarized or fragmented. Again, these 
results are not sensitive to the number of agents. 



3.4. Absorption time distribution 

The absorption time (number of steps until there are no more opinion interchanges) 
is of particular interest. With unbiased interactions, the expectation value of ab- 
sorption time can be determined using the results of Ross [23j as described in 
section 12.21 

Let ni,nj be the initial popularities of opinions i and j, respectively. The ex- 
pected number of interactions between these opinions is niUj. Then the expected 
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Fig. 4. Left: Cumulative distribution of absorption time for N = 98, Q = 7 and d = 1 . . . 6, under 
unbiased interaction. Points are simulation results and curves are nonlinear fits to the Generalized 
Extreme Value distribution. For all fits R? > 0.9996. Rigfit: Probability density for absorption 
time (derivative of the fitted curves on the left). 



E[tab.] = '^^^^^r '^ N\ (21) 



absorption time is given by the sum: 

E[tab.] = ^n,n,, (20) 

where the sum is over all compatible pairs of opinions such that i < j and 1 < 
\i-j\<d. 

When the initial distribution is uniform such that — N/Q for all i, this 
expression reduces to: 

d{2Q_d-l) 
2Q2 

Here, d{2Q — c? — l)/2 = c{Q,d) is the number of compatible opinion pairs for 
given Q and d. This number can be found recursively: It holds that c{Q,d) = 
c{Q, d — 1) + Q — d, because to go from tolerance d — 1 to tolerance d we need to 
connect opinion pairs (1, d), (2, d+1), . . . , (Q — d, Q), adding Q — d pairs to the sum. 
Using the initial condition c{Q, 1) = Q — 1 we get the result above. 

The expected absorption time scales as N'^, increases linearly in Q for fixed d, 
and increases monotonically with d. Narrow-minded agents (small d) require fewer 
steps to reach an absorbing state because they have fewer people to talk to. 

The absorption time distribution gives more information than the expectation 
value. Figure m shows the distribution with iV = 98 agents and Q = 7 opinions. The 
values are extracted from the same set of agent-based simulations used to produce 
Tabled 

With such wide distributions, it is more practical to use the Cumulative Dis- 
tribution Function (CDF) than to form a histogram (i.e., the probability density). 
Unlike a histogram, a CDF is always continuous, monotonically increasing, and not 
very sensitive to the number of runs used for averaging. Once a functional fit to the 
CDF is found, its derivative yields the probability density. 
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Table 2. Parameters for the GEV distribution determined with nonhnear fit. For most of the fits 
> 0.999. Bias is to pairwise majority. 



no bias 


a 


m 




1127 


1176 


0.2674 


1441 


2193 


0.1618 


454.1 


822.3 


0.1753 


811.1 


1466 


0.2287 


1197 


1966 


0.2156 


1483 


2455 


0.1407 


1593 


2834 


0.1035 


1610 


3040 


0.0960 



with bias 


a 


m 




46.03 


110.5 


-0.0096 


41.90 


164.6 


0.0735 


30.66 


118.2 


-0.0282 


37.56 


159.6 


0.0117 


41.65 


179.6 


0.0459 


46.32 


195.3 


0.0276 


48.00 


207.8 


0.0156 


46.38 


215.8 


0.0376 



The distribution of absorption times is very weh approximated by the General- 
ized Extreme Value (GEV) distribution |24j 1^ . The GEV distribution arises from 
the extreme value theory, a branch of statistics that deals with the maximum (or 
minimum) value from a set of i.i.d. random numbers. The theory is applied to model 
a wide variety of phenomena, including flood levels, rainfall, insurance claims, seis- 
mic events, and human life span. Since the absorbing state occurs for extremal 
popularity values (zero or maximum), it is reasonable that the GEV distribution 
appears in this problem, too. 

The GEV distribution has the cumulative distribution function: 



F(t; m, a, ^) = exp 



t — m 



(22) 



where the notation x+ indicates max(a;,0). The distribution has a location param- 
eter m, a scale parameter a > and a shape parameter ^. 

A minimum number of steps are required before any opinion is depleted of its 
followers; before that time the absorption probability is zero. Due to this lower limit 
on time, the GEV distribution slightly overestimates the simulation results at small 
time values {t < 1500). When biased interactions are used, the fit to GEV is almost 
perfect even at small time values (not shown in figure). 

Table [2] shows the parameters found by fitting the CDF of absorption time to the 
GEV distribution. The parameter is used to verify the goodness of fit; i?^ = 1 
indicates a perfect fit. The curve fits the simulation very well: The minimum value 
of R'^ is 0.9983 with Q = 3, d = 1, unbiased; for aU other cases R^ > 0.999. 



3.5. Survival probability of opinions 

Given a specific opinion q, one may ask how likely it is that this opinion will survive 
until the absorbing state. Table [3] shows the survival probability of opinions for 
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Table 3. Probability that a specific opinion survives in the absorbing state. Unbiased interactions. 



Q 


d 


3 


1 




2 


5 


1 




2 




3 




4 


7 


1 




2 




3 




4 




5 




6 



1 


2 


3 


4 


5 


6 


7 


0.6021 


0.3322 


0.6023 










0.3338 


0.3342 


0.3320 










0.5787 


0.3553 


0.4075 


0.3570 


0.5773 






0.4184 


0.2796 


0.2004 


0.2788 


0.4187 






0.3057 


0.2009 


0.2027 


0.1995 


0.3081 






0.2003 


0.2007 


0.2001 


0.1994 


0.1995 






0.5778 


0.3562 


0.3972 


0.3919 


0.3971 


0.3569 


0.5774 


0.4087 


0.2768 


0.2152 


0.2422 


0.2176 


0.2766 


0.4072 


0.3218 


0.2354 


0.1784 


0.1429 


0.1818 


0.2352 


0.3214 


0.2593 


0.1894 


0.1425 


0.1435 


0.1435 


0.1892 


0.2595 


0.2010 


0.1421 


0.1410 


0.1430 


0.1411 


0.1432 


0.2015 


0.1426 


0.1415 


0.1435 


0.1433 


0.1424 


0.1434 


0.1433 



(5 = 3,5 and 7, under unbiased interactions. The values are again obtained from the 
simulations used to produce Tables [T] and [2] Initial opinion distribution is uniform. 

Probabilities for a given Q and d do not add up to 1 because more than one 
opinion can be present in the absorbing state. For the cases where d = Q — 1 it 
can be predicted that the survival probability is 1/Q for each opinion [23] (see 
section [2?2|) . The values on the table are consistent with this prediction. 

The table shows that extreme opinions are more likely to be present in the 
absorbing state than central opinions. This should not be interpreted as agents being 
more likely to end up in extreme opinions. Actually, with unbiased interactions, the 
popularity of each opinion in the absorbing state is found to be approximately N/Q, 
independent of d or of the opinion's position on the spectrum (this does not hold 
when bias is present). 



3.6. Diffusion of opinions 

This section studies how fast an opinion spreads within a fully mixed population 
under unbiased interactions. 

Let nq{t) be the popularity of opinion q at time step t. By design, at each 
time step there is exactly one interaction between two compatible agents. If this 
interaction involves the opinion q, Uq will be changed by +1 or —1 with equal 
probability. If neither agent has opinion q, Uq is not changed. Therefore: 

t 

nq{t) = nq{0) + J2s^, (23) 

i=l 

where Si can be 1,0, or —1. Unlike the basic random walk, here the step size can be 
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zero. Define 

t 

SUq = Uq (t) ~nq{0)=Y^ S, . (24) 

The diffusion law, derived from the regular one-dimensional random walk, states 
that the variance of position is linear in time, where the constant of proportionality 
is twice the diffusion coefficient [24] . To arrive at a similar relation, consider the 
variance of Uq : 

(('^".)') = ( E E '^'^ ) = E E ' (25) 
/ i=i j=i 

where (• • • ) indicates an ensemble average (expectation value). Since Si and sj are 
independent random variables, their covariance (siSj) is {s^)Sij. Then all cross 
terms vanish, yielding 

(('5-,)')=E^?)=i^?)- (26) 

i=l 

The last step follows from the independence of successive steps. Therefore this 
process also obeys the diffusion equation with the diffusion coefficient D = (s^) /2. 

The next step is to calculate {sf), the expectation value of the change of pop- 
ularity Uq in one step. The value of is 1 if the interaction involves opinion q, 
and otherwise. Therefore (sf ) is the same as the probability that an interaction 
involves q. This value generally depends on q, on the number of opinions Q and on 
the confidence interval d. 

Define p{q; Q, d) to be the number of opinions that are compatible with a given 
opinion g. It is given by 

p{q; Q, d) = min(d, q - 1) + min(d, Q — q). (27) 

At every time step a compatible pair of opinions are selected for interaction. If q is 
one of these choices, there are p(g; Q, d) second choices. Therefore: 

{si) = ^^^'^'^''^ , (28) 

where the denominator is the number of all possible interactions, which is given in 
section [nH] as d{2Q — d — l)/2. With this substitution the diffusion coefficient is 
obtained as: 



^Q:Q:d = jtt^t; J • (29) 



min((i, q — 1) + mm{d, Q — q) 
d{2Q-d~ 1) 

For given Q and d, the diffusion coefficient increases linearly with the number 
of compatible opinions. Extreme opinions diffuse slowest because there are only d 
opinions compatible with them. In contrast, middle opinions diffuse fastest as there 
are 2d opinions compatible with them. 
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The result ((29)) indicates that, keeping q and Q the same, a larger value of the 
confidence interval d leads to faster diffusion for opinion q. If the number of opinions 
Q is increased while keeping q and d fixed, the diffusion of q slows down because 
there are more interactions which do not involve q. 

4. Identical agents on networks 

In fully mixed systems any agent can interact with any other agent directly. How- 
ever, such complete interaction networks are unrealistic for all but the simplest 
groups. In general, one can imagine each agent occupying a node in a graph and 
communicating with its neighbors only. Hierarchies can be represented by directed 
graphs where influence goes in one direction only. 

Scale- free networks or small- world networks are commonly used to simulate large 
communities [291 [1] [22]. Also, networks extracted from real- life relationships can be 
used [7]. 

A two-dimensional lattice would be an unrealistic representation of personal 
relationships. However, each agent can be taken to represent a geographic area. 
Then a lattice can be used to investigate spatial distribution of opinions and to 
simulate the spread (or decline) of a specific choice over a geographic region. The 
opinions may represent peaceful choices such as different brands of household items, 
or they may represent different armies engaged in hostilities. When two agents 
(areas) interact and one is converted to the opinion of the other, the brand or army 
associated with the winning opinion becomes dominant in both geographic areas. 

4.1. Absorbing states and blocking 

In the fully-mixed system, if the opinion q is a stationary opinion in an absorbing 
state, the popularities of all opinions compatible with q must be zero (Section E]). 
However, compatible opinions can coexist if the interaction network is not a com- 
plete graph. Two agents carrying compatible opinions may not be directly connected 
with each other, but to other agents with incompatible opinions. Figure [5] shows 
some examples of such static cases. Such arrangements would not be possible on 
fully connected interactions. The scarcity of connections isolates compatible opin- 
ions, creating pluralistic absorbing states. 

This blocking effect can be observed in a wide variety of interaction networks, as 
long as they are not complete graphs. Even if two agents are not directly connected, 
they can influence each other by means of other agents located between them. 
Therefore, an interaction becomes more likely and blocking less likely if, on average, 
there are only a few agents between any given pair. In other words, the strength 
of the blocking effect depends on the average distance between the nodes of the 
underlying graph [5T] . 

The coexistence of compatible opinions depends not only on the interaction 
network, but also on the confidence interval. Figure [5] illustrates how all opinions 
can coexist with Q — A, d — \. If, instead, the confidence bound d is taken to be 2, 
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(a) (b) (c) 




Fig. 5. Possible absorbing states with Q = 4 opinions and confidence interval d = I. (a) Static 
domains on a 2D lattice, labeled with opinions. Agents (lattice sites) with compatible opinions do 
not have a direct connection. For example, agents with opinion 1 and agents with opinion 2 are 
separated by agents with opinion 4. (b) and (c) show some static states on different interaction 
networks. The circles are the agents and the labels indicate opinions. 

such perfect pluralism is impossible, regardless of the connections. This is because 
opinions 2 and 3 are compatible with all opinions, and their interactions cannot be 
blocked. 

In general, if Q — d < 1 + d holds, then opinions labeled Q — d through 1 + d 

are connected to every opinion. In that case, any absorbing state can have at most 
2(Q — d — 1) stationary opinions (this is not the opinion value, but the count of 
opinions) . 

5. Conclusions 

This study presents a new bounded-confidence model for discrete opinions. The 
model addresses situations where a compromise is not possible; in any interaction 
one of the agents, selected randomly, converts the other to its opinion. Similar 
models exist, but as far as we know, this is the only model that combines arbitrary 
number of discrete opinions, bounded confidence, and lack of compromise. 

If the population is fully-mixed and if agents are eqtially likely to convert each 
other, the model is equivalent to the multi-player gambler's ruin problem. Using 
this equivalence it is shown that the system always ends up in an absorbing state. 
The time (number of opinion exchanges) required to reach that state scales as N"^ 
and increases monotonically with the confidence bound d. 

An absorbing state consists of at least one, at most [{Q — 1) / {d+l)}+l surviving 
opinions. The probability of consensus, for a given Q, decreases with the confidence 
bound d. 

In fully-mixed systems, the distribution of absorption time is found to agree 
very well with the Generalized Extreme Value distribution. 

In fully mixed systems with imbiased interactions, the spread of the popularity 
of an opinion obeys the difFtision equation. The difFtision coefficient of an opinion is 
proportional to the number of opinions compatible with that opinion. 

When the interaction network is not a complete graph, a given opinion may 
become stationary even if compatible opinions survive. This is possible if the path 
between compatible agents is blocked by incompatible agents. 
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The model is flexible and it can be extended in different ways. The agent pop- 
ulation can be made heterogeneous with varying tolerances, mixing tolerant and 
stubborn agents. Agents may carry vector opinions and may interact with each 
other only if they are compatible on all dimensions. Noise can be introduced by 
making agents switch to another opinion spontaneously. Interaction bias can be 
based on individual properties of agents, such as opinion values. Agents can keep a 
record of their interaction history and/or opinions of agents they meet (compatible 
or not) so that the conversion probability can be based on that history. 

On the macroscopic level, the interaction network of agents may be changed, 
and properties of diffusion and absorbing states can be compared with the fully- 
mixed case. Some obvious modifications include higher dimensional lattices, random 
networks, small-world networks or scale-free networks. In such networks, the distri- 
bution of opinion clusters can be investigated. 
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